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[57] ABSTRACT 

A technique for reducing the spurious signal content in 
digital sinusoid synthesis. Spur reduction is accomplished 
through dithering both amplitude and phase values prior to 
word-length reduction. The analytical approach developed 
for analog quantization is used to produce new bounds on 
spur performance in these dithered systems. Amplitude 
dithering allows output word-length reduction without intro- 
ducing additional spurs. Effects of periodic dither similar to 
that produced by pseudo-noise (PN) generator are analyzed. 
This phase dithering method provides a spur reduction of 
6(M+1) dB per phase bit when the dither consists of M 
uniform variates. While the spur reduction is at the expense 
of an increase in system noise, the noise power can be made 
white, making the power spectral density small. This tech- 
nique permits the use of a smaller number of phase bits 
addressing sinusoid look-up tables, resulting in an exponen- 
tial decrease in system complexity. Amplitude dithering 
allows the use of less complicated multipliers and narrower 
data paths in purely digital applications, as well as the use 
of coarse-resolution, highly-linear digital-to- analog convert- 
ers (DACs) to obtain spur performance limited by the DAC 
linearity rather than its resolution. 

10 Claims, 7 Drawing Sheets 
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METHOD AND APPARATUS FOR 
SPUR-REDUCED DIGITAL SINUSOID 
SYNTHESIS 


ORIGIN OF INVENTION 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the io 
Contractor has elected not to retain title. 

TECHNICAL FIELD 

The present invention relates to a method and apparatus 15 
for adding a random noise sequence to the phase and/or 
amplitude representation of a digitally generated sinusoid to 
reduce the required number of bits without greatly increas- 
ing the spurious content. The invention is useful for gener- 
ating spectrally pure digital sinusoids without employing a 2 o 
large number of bits in phase representation, which would 
otherwise result in exponentially large look-up tables, or, in 
amplitude representation which, in turn, would result in 
increasing the complexity of multiplications and in limiting 
the bandwidth of digital-to-analog conversion. The inven- 2 s 
tion adds a random noise sequence to the phase and/or 
amplitude samples as appropriate, followed by a rounding or 
reduction of the number of bits originally in the represen- 
tation. This addition of noise “dithers” the resulting phase 
and/or amplitude values, reducing the spurious signal con- 
tent of the smaller wordlength representation. 

BACKGROUND ART 

It is well-known that adding a dither signal to a desired 
signal prior to quantization can render the quantizer error 
independent of the desired signal. Classic examples of this 
work deal with the quantization of analog signals. Advances 
in digital signal processing speed and large scale integration 
have led to the development of all-digital receiver systems, 4Q 
direct digital frequency synthesizers and direct digital arbi- 
trary waveform synthesizers. In all these applications, 
because finite word-length effects are a major factor in 
system complexity, they may ultimately determine whether 
it is efficient to digitally implement a system with a particu- 45 
lar set of specifications. Earlier work has presented a tech- 
nique for reducing the complexity of digital oscillators 
through phase dithering with the claim of increased fre- 
quency resolution. Recent research has suggested mitigation 
of finite- word-length effects in the synthesis of oversampled 50 
sinusoids through noise shaping. It would be useful if the 
analysis techniques used for quantization of analog signals 
can be applied to overcome finite-word-length effects in 
digital systems. It would also be advantageous if appropriate 
dither signals can be used to reduce word lengths in digital 55 
sinusoid synthesis without suffering the normal penalties in 
spurious signal performance. 

Conventional methods of digital sinusoid generation, e.g., 
FIG. 1, result in spurious harmonics (spurs) which are 
caused by finite word-length representations of both ampli- 60 
tude and phase samples. Because both the phase and ampli- 
tude samples are periodic sequences, their finite word-length 
representations contain periodic error sequences, which 
cause spurs. The spur signal levels are approximately 6 dB 
per bit of representation below the desired sinusoidal signal. 65 

A search of the most relevant prior art resulted in the 
following U.S. Patents: 


4,652,832 

Jasper 

4,926,130 

Weaver 

4,994,803 

Blackham 

5,014,231 

Reinhardt et al 

5,017,880 

Dugan et al 

5,029,120 

Brodeur et al 

5,036,294 

McCaslin 

5,073,869 

Bjerede 

5,091,921 

Minami 

5,108,182 

Murphy 

5,121,409 

Goss 

5,128,623 

Gilmore 

5,131,750 

Gravel et al 

5,151,661 

Caldwell et al 

5,162,746 

Ghoshal 

5,166,629 

Watkins et al 

5,184,093 

Itoh et al 


Of the foregoing patents, the following appear to be the 
most pertinent: 

U.S. Pat. No. 5,073,869 to Bjerede is directed to suppres- 
sion of spurious frequency components in a direct digital 
frequency synthesizer. The improved direct digital fre- 
quency synthesizer has a coarse 48/fine accumulator 40 and 
a non-linear digital-to-analog converter 42. The accumulator 
suppresses the generation of spurious frequency components 
in the analog waveform output of the synthesizer by ran- 
domly dithering the phase at which the coarse component 
accumulator 48 is incremented by the phase accumulator 40. 
The sample-and-hold circuit 44 mitigates spurs induced by 
30 glitches from the digital to analog converter by sampling the 
analog signal 56 at times when the glitches are not present. 

U.S. Pat. No. 4,994,803 to Blackham is directed to a 
random number dither circuit for digital to analog output 
signal linearity. Distortion is reduced in a linear circuit by 
adding a random digital number to the input and subtracting 
the equivalent random analog number from the output. A 
digital number from random number generator 14 is 
summed with the input digital number and is converted to an 
analog number by DAC 18. The summed input and random 
number is converted to an analog signal by DAC 24 and the 
output of 18S is subtracted in summer 22 to provide the 
distortion-free output. 

U.S. Pat. No. 5,162,746 to Ghoshal is directed to a circuit 
for attenuating phase jitter in a clock signal that includes a 
dithering circuit, a phase locked loop, and a digitally con- 
trolled oscillator. The dithering circuit modulates the digi- 
tally controlled oscillator to improve rejection behavior 
when the incoming clock frequency is substantially the same 
as the divided down oscillator signal. 

U.S. Pat. No. 5,108,182 to Murphy is directed to a digital 
path length control for a ring laser gyro comprising a dither 
counter, a dither DAC, a control counter, a control DAC, and 
a piezoelectric transducer path length control driven by two 
DAC’s. The dither modulates the operating conditions to 
center rather than at either extreme of the system error. 

U.S. Pat. No. 5,036,294 McCaslin is directed to a phase 
locked loop having low frequency jitter compensation. An 
output clock that is in-phase with a reference clock utilizes 
a dither circuit to control the switching of the phase locked 
loop. A switch capacitor phase locked loop 10 has a dither 
portion 16 that translates the low frequency intrinsic jitter to 
a higher frequency jitter where it is more tolerated. 

None of the aforementioned prior art discloses a method 
and apparatus for exponentially decreasing the complexity 
of numerically-controlled oscillators and direct digital fre- 
quency synthesis with only a small increase in system noise. 
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STATEMENT OF THE INVENTION 

The technique of the present invention reduces the rep- 
resentation word length without increasing spur magnitudes 
by first adding a low-level random noise, or dither, signal to 5 
the amplitude and/or the phase samples, which are originally 
expressed in a longer word length. The resulting sum, a 
dithered phase or amplitude value, is truncated or rounded to 
the smaller, desired word length. Of course, either the 
amplitude or the phase or both can be dithered. In phase 1Q 
dithering the spurious response is determined by the type of 
dithering signal employed. In amplitude dithering the spu- 
rious response is determined by the original, longer word 
length. While the amplitude-related spurious is generally 
related to the phase-related spurious, the pre-dither ampli- ]5 
tude word length is made long enough to satisfy spur power 
specifications. Then the exact relationship is unimportant, 
and since the phase dither signal is independent of the 
amplitude dither signal, the amplitude and phase dithering 
processes can be treated independently. 20 

The applicant first describes the quantizer model. Ampli- 
tude and phase quantization effects are then reviewed, and 
simple new bounds on spurious performance are presented. 

In contrast to bounds in the prior art, the new bounds are 
straight-forward and require little information about the 25 
signal to be quantized. The derivations of the new bounds 
provide motivation for new analysis of dithered quantizer 
performance. An analysis of dithering with a periodic noise 
source is presented. The periodic noise source is considered 
because of its similarity to implementations involving linear 30 
feedback shift registers (LFSRs), or Pseudo-Noise (PN) 
generators. New analysis of phase dithering effects is pre- 
sented, followed by simulation results and a design example. 

OBJECTS OF THE INVENTION 35 

It is therefore a principal object of the invention to provide 
a method and apparatus for implementing high performance 
digital synthesizers with dramatically reduced complexity 
and while avoiding high noise power spectral density levels. 

40 

It is another object of the invention to provide a method 
and apparatus for using a digital dithering approach to spur 
reduction and reduction in the complexity of numerically- 
controlled oscillators, frequency synthesizers and other peri- 
odic waveform generators. ^ 

It is still an additional object of the present invention to 
employ a new approach to amplitude and/or phase dithering 
which permits the reduction of word length while reducing 
spurs thereby facilitating significant reductions in the com- 
plexity of numerically-controlled oscillators, frequency syn- 50 
thesizers and other periodic waveform generators. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The aforementioned objects and advantages of the present 
invention, as well as additional objects and advantages 55 
thereof, will be more fully understood hereinafter as a result 
of a detailed description of a preferred embodiment when 
taken in conjunction with the following drawings in which: 
FIG. 1 is a block diagram of a conventional digital 
sinusoid generator; 60 

FIG. 2 is a graphical representation of the input/output 
relation of a midtread quantizer; 

FIG. 3 is a conceptual block diagram of a waveform 
generator; 65 

FIG. 4 is a block diagram of a two-step waveform 
generator; 


FIG. 5 is a block diagram of a uniform dithered quantizer; 

FIG. 6 is a block diagram of a system for 1 8 dBc per phase 
bit spur reduction; 

FIG. 7 is a graphical representation of a power spectrum 
of an 8 sample/cycle sine wave without dithering; 

FIG. 8 is a graphical representation of a power spectrum 
of an 8 sample/cycle sine with amplitude dithering; 

FIG. 9 is a graphical representation of a power spectrum 
of a 5-bit phase- truncated sine wave without phase dither- 
ing; 

FIG. 10 is a graphical representation of a power spectrum 
of a 5-bit phase-truncated sine wave with first-order phase 
dithering; 

FIG. 11 is a graphical representation of a power spectrum 
of a 5-bit phase-truncated sine wave with second-order 
phase dithering; 

FIG. 12 is a graphical representation of a worst-case 
power spectrum of a sinusoid with first-order phase dither- 
ing and amplitude dithering; and 

FIG. 13 is a block diagram of a spur-reduced direct digital 
frequency synthesizer in accordance with an embodiment of 
the present invention. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT 


Quantizer Model 

When a discrete-time input signal x[n], is passed through 
a uniform mid-tread quantizer, the output signal, y[n], can 
always be expressed as y[n]-x[n]+e[n] where e[n] is the 
quantization error, a deterministic function of x [n]. The 
input to the quantizer is mapped to one of the 2 b levels, 
where b is the number of bits which digitally represent the 
input sample. Output levels are separated by one quantizer 
step size, A=2~ b . Throughout this description A a will be used 
as the step size for amplitude quantization results, A P will be 
used for phase quantization results, and A will be used if the 
result applies to both amplitude and phase quantization. 
Similar subscripting will be used on the quantization error. 

The input/output relation of a mid-tread quantizer appears 
in FIG. 2. If the input does not saturate the quantizer then the 
quantizer error is: 





j2nkx[n ] 
A 


( 1 ) 


If the input signal is bounded so that lx[n]I^A Q where A Q = 
V^-A, then the quantizer does not saturate and le[n]I^A/2. 
Throughout this paper, quantizers are always operating in 
non-saturation mode. 

Amplitude Quantization Effects 

Let a discrete-time sinusoid with amplitude A^A e and 
frequency co 0 be the input to a mid-tread quantizer. If the 
sinusoid is generated in a synchronous discrete-time system, 
co 0 can be expressed as 2n times the ratio of two integers. 
The input sequence is the periodic with a finite period. Since 
the error sequence, e A [n], is a deterministic function of the 
input sequence, it is periodic with a finite period as well. So 
e A [n] has a Fourier series representation as a function of 
time. Therefore, the spectrum of the error sequence will 
consist of discrete frequency components (spurs) which 
contaminate the spectrum of x[n]. 

The following argument leads to an upper bound on the 
size of the largest frequency component in the spectrum of 
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e A [n]. Assuming the quantizer is not saturated by the input 
signal x[n], the maximum possible quantization error is 
A a / 2, where A a is the amplitude quantization step size. The 
total power in e A [n] is then bounded by A a 2 /4. By ParsevaF s 
relation, the sum of the spur powers in the spectrum of e A [n] 5 
equals the power in e A [n]. In order to maximize the power 
in a given spur the total number of spurs must be minimized. 
Since e A [n] is real with the exclusion of static (DC) spurs 
and spurs at half the sampling rate, DC offsets and half 
sampling rate spurs can be corrected by appropriate calibra- 10 
tion and filtering, the maximum power in a spur occurs when 
there are two frequency components at -ko 3 and -ft) 3 , with 
equal power. With two frequency components, e A [n] is 
sinusoidal. Therefore since the power in e A [n] is ^ A a 2 /4, the 
power in a single spur is £§A a 2 /8. 15 

Since x[n] is real, its spectrum consists of two frequency 
components, at +co 0 and -co 0 , each having power A 2 /4. 
Using the above bound on spur power, the spur to signal 
ratio (SpSR) is ^A a 2 /(2A 2 ). If A=A Q ~ l /2 provided b is not 
small, then in decibels with respect to the carrier (dBc), 20 
SpSR ^3-6 b dBc, where A a = 2~*, and b is the word length 
in bits. In summary, the upper bound above on power in a 
spur caused by amplitude quantization exhibits -6 dBc per 
bit behavior. 

Phase Quantization Effects 25 

Now let a phase waveform, cj>[n], be the input to the 
mid-tread quantizer. The phase waveform, <j>[n]=(fn+<j)/27t) is 
a sampled sawtooth with amplitude ranging from 0 to 1. The 
fractional operator, (x), is defined so that (x)=x mod 1, e.g., 
(1.3)=0.3. Since <J>[n] is generated by a synchronous, finite- 30 
word-length, discrete-time system, it has a finite period. The 
signal output from the quantizer can be expressed as <j>[n]+ 
ejn], where e p [n] represents the error introduced by quan- 
tization. Again, since <j>[n] is periodic, e p [n] is periodic with 
a period less than or equal to the period of <)>[n]. After 35 
multiplication by 2 k and passage through the ideal function 
generator, the output signal is y[n]=A cos(27c<(> [n]+2rce p [n]). 

If the quantizer has many levels, i.e., >16, then e p [n]«l. 

The approximation y[n]~A cos (2 mj> [n])-2nAe p [n]sin(2rc<|) 

[n]) is obtained using small angle approximation for cosine. 40 

Since ejn] and <j>[n] are periodic, the total error 2A7te p 
[n]sin(2rc<j> [n]) is periodic. The total error is a real signal, so 
using the above amplitude quantization argument, the maxi- 
mum spur power is bounded by (max{27tAe p [n]sin(2n;<j>[n] 
)} 2 )/2. This equals 7t 2 A 2 A p 2 /2, where A p =2~ b and b bits are 45 
used to represent phase samples. By the above approxima- 
tion for yfn] and the bound on the spur power SpSr^2it 2 
A p 2 =13-6 b dBc, independent of the signal amplitude, A. In 
practice this simple bound can be improved by about 9 dB, 
but it demonstrates the maximum spur power behavior of -6 50 
dBc per phase bit. 

Amplitude Dithering 

It will now be shown that rounding the sum of an already 
quantized sinusoid and using an appropriate dither signal 
cause spurious magnitudes which depend on the original 55 
(longer) word length, not the output (shorter) word length. 
This phenomenon occurs at the expense of increased system 
noise from the addition of the dithering signal. An important 
finite word-length dithering system is subsequently shown 
to be equivalent to the continuous-amplitude uniformly- 60 
dithered system. 

Consider the conceptual block diagram for a waveform 
generator shown in HG. 3. The b-bit quantizer can be split 
into two parts: a high resolution B-bit quantizer (B>b) 
followed by a unit for truncation or round to b bits. The 65 
resulting waveform generation system is shown in FIG. 4. 
Thus, the generation process consists of two separate steps: 


production of a high-resolution waveform and reduction of 
the word length. The number of bits used to represent the 
high-resolution samples should be sufficient to guarantee the 
desired spectral purity. Then, the word length should be 
reduced without creating excess signal-dependent quantiza- 
tion error. 

The input in FIG. 5 xfn], is a B-bit representation of A 
sin(27t<j>[n]), where <))[n] is the phase value at time n. The 
input can be expressed as x[n]-Asin(27C<j>[n])+e AC) [n3, where 
e AO [n] is the quantization error due to some previous quan- 
tization of A sin(2n<|>[n]). For example, x [n] could be the 
output from a sine look-up table with output word length B 
bits. The dither signal, zjn], is white noise uniformly 
distributed in [-A a /2, A a / 2), where A A =2~ b . 

After zjn] and xfn] are added, the output is rounded to 
retain only the b most significant bits. The rounding can be 
modeled as a uniform quantizer with step size A a . The 
amplitude A is chosen to avoid saturating this quantizer 
when the dither signal is added, i.e., A+A a /2^A q . 

The output from the quantizer can be expressed as y[n] 
=x[n]+zJn]-fe A [n]. The characteristic function of the dither 
signal, z M [n], is: 


F z ( a) = £{exp(/az[rt])} = sine 


aA a 
2k 


(2) 


where sine (x)=sin(7tx)/7tx. Since the characteristic function 
of zjn] has zeros at integer multiples of 27t/A A , e A [n] will be 
white and uniformly distributed over [~A a /2, A a / 2). There- 
fore the quantization error e A [n] is a white, wide- sense 
stationary process that does not contribute spurious harmon- 
ics to the output spectrum of y[n]. 

The only spurious components in y[n] are due to e AO [n], 
which contains the spurious components originally in xfn]. 
It remains to comment on the noise, i.e., the power not 
isolated in discrete spurious frequency components. If the 
sequences e A [n] and zjn] are uncorrelated, adding the 
variances of the quantization error, A a 2 /12, and the dither 
process, A a 2 /12 yields a white noise power of A a 2 /6. This 
approximation is twice the variance of a quantization system 
with no dithering signal. Note that e AO [n] also contributes a 
white noise term, and that, in general, e A [n] and zjn] are not 
uncorrelated. However, these two effects are dominated by 
the A a 2 /6 behavior of the white noise power. In summary, 
y[n], which is quantized to b bits, exhibits spurious perfor- 
mance as if it was quantized to B bits (B>b), at the expense 
of doubling the white noise power. 

Because the input xfn] is expressed as a B-bit value, an 
important equivalent system to continuous-amplitude, uni- 
formly-dithered word-length reduction can be constructed 
by replacing the uniformly distributed dither signal, zjn], 
by a finite word-length of B-b bits, and is said to be 
discretely and evenly distributed over the quantized values 
in the region f-A A /2, A a / 2). Heuristically, zfn] randomizes 
the portion of the finite word-length input, x [n], that is about 
to be thrown away by the rounded truncation. This process 
is equivalent to continuous uniform dithering, since if xfn] 
is padded out to an infinite number of bits by placing zeros 
beyond the least significant bit (LSb), then only the B-b 
most significant bits of zjn] will have an effect on the 
resulting sum x[n]+zjn]. All of the bits below the most 
significant B-b are added to zero, and cannot beget a carry. 
The output, yfn], is identical in both systems. Therefore 
zjn], continuously, uniformly distributed over [-A a /2, 
A a / 2) can be replaced by zfn], and yield the same spurious 
response for yfn]. 

It appears that the finite word-length dither signal, zfn], 
could be generated by a linear feedback shift register 
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(LFSR), or PN generator. This will be strictly true if and 
only if the PN generator has an infinite period, since at this 
time, the dither signal is required to be white. However, it is 
not surprising that the ideal behavior is approached as the 
period of the PN generator gets longer. With a sufficiently 5 
long period, the case where spur magnitudes are limited by 
the original word length can be achieved. A simple model for 
a system implementation using a periodic random sequence 
which can be approximated by a PN generator will now be jo 
described. 

Effect of Periodic Dither 

The use of a periodic dither signal with a long period, L, 
for both amplitude and phase dithering is now analyzed. 
Since the dither signal is periodic, the discrete frequency 15 
components in its spectrum will contaminate the desired 
signal. It is shown that the period can be chosen to satisfy 
worst case spurious specifications. The case where the dither 
signal is generated using one uniform variate (M=l) is 2 o 
given. When the dither signal is the sum of M independent 
uniform variates (M>1), the analysis is the same because the 
resulting signal is an i.i.d. sequence of random variables. 

Instead of using the white dither process, zjn], described 
in the previous section, consider a substitute, zjn]. The 25 
dither process, zjn] is periodic with period L. Any two 
samples, zjn] and zjn+m], where m*0 rood L, are inde- 
pendent. Samples of zjn] are uniformly distributed between 
[-A/2, A/2), and the quantization step size is A. 30 

When zjn] is used as the dither signal, let the quantizer 
error be called ejn]. The autocorrelation of zjn] when the 
lag m, is an integer multiple of Lis equal to R 2Z ^ L [0]=A 2 /12. 

In the PN generator approximation to this noise source, 
L=2 rf -1 where d is the length of the shift register in bits. At. 35 
other lag values, the samples of zjn] are independent, and 
since they have zero mean, the autocorrelation is zero. 
Therefore: 

40 

fi W m] = ~fr 5[m mod L] = ^ -mr exp ( j 2 t }L ) 

Therefore, the spectrum of zjn] consists of L discrete 45 
frequency components with power A 2 /12L. 

In the autocorrelation expression for ejn], the expecta- 
tion is taken over the random variables zjn] and zjn+m]: 


R e L e L [n,n + m] = 


50 


( 3 ) 


Z X a*[n](Xi*[n + 

k=-oo l——oa 

teQ teQ 

55 

m]E | exp | (kzdn] ~ kdn + m]) J J 


where: 


a kin] = 


A (- 1 ? 
j2nk 


exp 


j2nks[n] 

A 


60 


The desired signal to which the dither signal zjn] is added 
is s[n]. Using the notation from earlier sections, in phase 65 
quantization, s[n]=<))[n] and in amplitude quantization s[n]= 
x[n]. When the lag is not an integer multiple of L, 


E | exp | (kzdn] - Izdn + m]) j j= 


f C .. D ( \ 

L f -flrtzrtn + m] 


\ exp l A 1 

\ E \ <*V\ A 

I 


Fz(-2f- )J — X— )=8M8[Q. 

This last fact is true because the characteristic function of 
zjn] has zeros at all non-zero integer multiples of 2nJA 
(Equation 2). But since the sums over k and 1 never assume 
the value 0, the autocorrelation function is zero when the lag 
is not 0 mod L. When the lag is 0 mod L: 

E { exp | (kzdn] - Izdn + m]) j j= 

„ \ I j2n(k~l)zdn] ) \ wpi n 

E | exp ^ j J j=8[*-J]. 

This results in: 

R e L e L [n,n + m] = 


A 2 " 1 ( 2nk , , , , ^ n \ 

S- S — 9— cos — r — (s[n] + 

2k k= 1 k \ ** / 

Setting m=0 in (4) yields the power in ejn]: 


r Vl m = 



“ A 2 

h=\ k 2 12 


From Equation 4, ejn] is a cyclo- stationary process because 
s[n] has a finite period, N. Using the results of Ljung, 
spectral information is obtained when Equation 4 is aver- 
aged over time. Note that when the lag, m, is not only an 
integer multiple of L, the period of the either, but also an 
integer multiple of N, the autocorrelation function equals 
A 2 /12, independent of n. The smallest nonzero lag that 
satisfies these two conditions is the least common multiple 
of L and N, denoted by qL where q is an integer. Therefore, 
the period of the time-averaged author relation function, 
R«Jm]=Avg (R eZ ^Jn, n+ml), is at least L aria at most qL. 
Let the period equal cL, where c is an integer, l^c^q. The 
function R e Jm] can be expressed as a sum of cL weighted 
complex exponentials: 


R e L [m ] = p\ exp | Jl—rL. j, m = . . . , -1,0,1, 2, . . . 

where 


P 1 = 


cL-1 _ 

£ R e r [m]exp 
m = 0 L 


j2nml 

cL 




_L_ 

cL 


T R e [nL\Qx p 
n=0 L 


j2nnl 

c 


The last equality is true since the autocorrelation function in 
Equation 3 and its time-average, R e Jm] are zero for lags not 
equal to integer multiples of L. The weights, P l5 are the 
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power magnitudes of the spurs. Since R e Jm] is less than or 
equal to A 2 / 12, the spur power can be bounded: 

< A 2 < A 2 

Pl = YlcL = 12L * 


10 

and <j)[n] is periodic. At all sample times n, E{€[n]}=-A p 
p[n](l-p[n])+A p (l“p[n])p[n]=0. 

Information about the spurs and noise in the power 
spectrum of x[n] is obtained from the autocorrelation func- 
5 tion. The autocorrelation of x[n] is: 


Equality is achieved when the period of the time-averaged 
autocorrelation function is exactly L, the period of the dither. 

As L-xso, the spacing between spurs goes to zero in the 10 
spectra of both ejn] and zjn]. The power in an individual 
spur goes to zero, but the density (power per unit of 
frequency) tends to a constant. Thus, ideal white noise 
behavior is approached. While z L [n] and ejn] are correlated 
in general, the worst case spur power scenario coherently 15 
adds the power spectra from both processes. For this reason, 

L should be chosen to satisfy 2A 2 /12L<P mac , where P max is 
the maximum acceptable spur power. When constructing a 
dither signal as the sum of M^l independent, uniform 
variates the noise autocorrelation becomes R z£zL [m]=MA 2 / 20 
12S[m mod L]. The analysis follows closely to that for M=1 , 
and L should be chosen to satisfy (M+l)A 2 /12L<P max . 

Since the desired signal has finite word length, it is 
equivalent to round or truncate the dither signal to an 
appropriate finite word length. The truncated periodic noise 25 
source is an approximation to the behavior of an implemen- 
tation using a PN generator which produces a periodic 
sequence of discretely and evenly distributed random num- 
bers. 

Phase Dithering 30 

Phase dithering is now analyzed using a continuous, 
zero-mean, wide-sense stationary sequence. As described 
above in reference to amplitude dithering, an evenly distrib- 
uted discrete random sequence is equivalent to continuous 
uniform dithering when the initial phase word is quantized 35 
to a finite number of bits. 

Let the digital sinusoid to be generated be: 

T[n]=cos(2rc$[n]) (5) 

and let the phase process be given by: 

0[n]=4[n]+£[n]=tfn-f4>/(27i))+€[n] (6) 

so that the desired phase is <j>[n], measured in cycles with a 
frequency of f cycles per sample, and a static offset of 
radians. The total quantization noise is e[n]=e p [n]+z[n], the 
sum of the dither signal and the quantizer error. Using small 
angle approximations: 

T[nl^os(27cfn+0)-2jce[n]sin(27cfn+O)+O((max(£[nl)) 2 ). 

The total quantization noise will be examined by consider- 
ing the first two terms above, and then the second-order, 
0([max(e[n])] 2 ), effect. 

A) FIRST ORDER ANALYSIS 

Since the quantization error after dithering is independent 
of the input signal e[n] is uncorrelated with the desired 
sinusoids. Without loss of generality, and for ease of nota- 
tion, let us shift the uniformly distributed dither random 
variate range to (0,A p ). The total phase quantization noise 
e[n] will be e[n]=-A p p[n] with probability (l-p[n]), and 
€[n]=(l-p[n])A / , with probability p[n]. The value p[n] is the 
distance from the initial high-precision phase value, <j>[n], to 
the nearest greater quantized value normalized by the phase 
quantization step size A p . The value of the probability 
sequence p[n] varies periodically, since p[n]=<j>[n] mod A p , 


£{jt(n]*[n + m]}~ cos(2nJn + O)cos(2ic/(n + m) + <t») + 

47i 2 sin(27Efn + 0)sin(2 nj{n + w) + 

0)E{€[n]€[n + m]} + 0(Ap4). 

In the absence of error (e[n] identically zero), the autocor- 
relation of the process x[n] is 

E{'c(n]'i[n+m3}=cos(2TtfiH-0)cos(2Ttf[n+mH<5) 

which is equal to the first term in the autocorrelation when 
the error sequence is non-zero. The second term in the 
autocorrelation is due entirely to the error. 

When e[n] is not identically zero spectral information is 
obtained by averaging over time (9), resulting in: 

R rt [m]-( 1 /2)[l+4jc 2 R i; Jm]]cos( 2itfm) 

where R ee [m]=Avg n (E{€[n]€[n+m]}), time-averaged auto- 
correlation of the total quantization noise. 

The power spectrum of x[n], the Fourier transform of the 
autocorrelation, is the power spectrum of desired sinusoid of 
frequency f with the total quantization noise amplitude 
modulated on the desired sinusoid. Note that since R ee [m] 
=0(A p 2 ), and A p «l, the modulation index is small. 

To a first-order the AM signal produced by phase dither- 
ing is clear of spurious harmonics down to the level due to 
periodicities in the dither sequence. The next section will 
examine spur performance in more detail, but first it is 
important to consider the noise power spectral density added 
in the phase dithering process. This result is achieved by 
deriving the signal-to-noise ratio (SNR). 

Recall that for any fixed time n, the probability distribu- 
tion of e[n], a function of p[n], is determined by the input 
signal, but the outcome of e[n] is determined entirely by the 
outcome of the dither signal z[n]. When z[n] and z[n+m] are 
independent random variables for non-zero lag m, efn] and 
e[n+m] are also independent for im*0, and hence e[n] is 
spectrally white. In this case, the autocorrelation becomes: 

^[mj^VzcosClTtfmyt-ZTriSfmjVarCe) 

where Var(e) is the time-averaged variance of the total 
quantization noise, and 8[m] is the Kronecker delta function 
($[0]=1, S[m]=0, m^O). The resulting SNR is: 

SNR-( 1 /2)/(27c 2 Var(e))-l/(47t 2 Var(e)). 

When the dither signal is constructed from one uniform 
[-Ap/2, Ap/2) random variate, the error e[n] is bounded 
between ~A p to A p with A p =2~ b , and b is the number of bits 
in the phase representation after the word-length is reduced. 
The number of bits, b, must be large enough to satisfy the 
small angle assumption earlier, typically bi?4. The time- 
averaged variance of e[n] is less than or equal to 2~ 2b I4, and 
the SNR is: 

SNR 1 2 2b * 4/(4tc 2 ) = (202?)log lo (2) - 201og lo (7t)dB § 

6.022? - 9.94dB. 

Since the sinusoid generated is a real signal, the signal 
power in the SNR will be divided between positive and 
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negative frequency components. If the sinusoid is the result 
of a discrete-time random process with sampling frequency 
fy, then the resulting noise power spectral density (NPSD) 
will be given by: 

5 

NPSD = —[SNR/2 + 10 log 10 (£/2)]dBc/Hz ^ (7) 

6.93 - 10 log 10 (£/2) - 6.02 b dBc/Hz. 

As an example, TABLE I gives noise power spectral Q 
densities as a function of the number of bits per cycle, b, at 
a 160 MHz sampling rate, calculated according to the above 
formula. 


TABLE I 


b (bits/cycle) 

Noise Power Spectral Density 

5 

-102.20 dBc/Hz 

6 

-108.22 dBc/Hz 

7 

-114.24 dBc/Hz 

8 

-120.26 dBc/Hz 

9 

-126.28 dBc/Hz 

10 

-132.30 dBc/Hz 

11 

-138.32 dBc/Hz 

12 

-144.35 dBc/Hz 


25 

B) SECOND ORDER ANALYSIS: RESIDUAL 
SPURS 

For a worst-case analysis of second order effects, expand 
the initial cosine from Equation 5 and 6 by the sum of angles 30 
formula: 

T[n]=cos(2Tce[n])cos(27cj 7 n-Kl>)-sin(27iefn])sin( 2?rf n+cD). 

Information about the spurs in the power spectrum of x[n] is 35 
obtained from the autocorrelation function at nonzero lags. 
When the dither sequence, z[n], is sequence of i.i.d. variates, 
the autocorrelation function for x[n], with lag m not equal to 
zero, is: 

R rt [n,n+m]=E{x[n] x[n-Hn] }=E{ x[n] }E{x[n+m] }. 40 

The expected value of x[n] is a deterministic function of 
time. From the above expression, it follows that spectral 
information about the random process x[n], with the excep- 
tion of noise floor information, is contained in E{x[n]}. Call 45 
the sequenced of E{x[n]} the “expected waveform”: 

E{x[m ] } = £{cos(27ie[ro])}cos(27t/m + O) - 

£{sin(27te[ro])}sin(2itym + O) = ^ 

(1 - 2jr 2 E{e 2 [m]})cos(27i/m + O) - 

27t£{e[m]}sin(27t/m + 3>) + 0(Ap 3 )- 

Since e is zero mean at all sample times this reduces to: 55 

E{T:[m]}=(l-2jt 2 E{6 2 [m]})cos(27cfm+0)+0(A/). 

It remains to consider E{e 2 [m]}, which we evaluate by 
using the probability sequence p[m] from the previous 6Q 
section: 

E{e 2 [m] }=A p 2 p 2 [m](l-p[m])+Ap 2 (l-p[m]) 2 p[m]=A/(p[m3-p 2 [m] 

)• 

Since p[m] is bounded between 0 and 1, the function 65 
u[m]=p[m]-p 2 [m] is bounded between 0 and l A, with its 
maximum value of Va at pfmj^A 


12 

Since u[m] is bounded between 0 and l A, it must have 
some non-zero DC (average) component. Any remaining 
components can be periodic in the worst case. Since all 
non-linear operations have been performed, conservation of 
power (energy) arguments can be used to determine the total 
non-DC error power. The total power in the DC component 
of u[m] is equal to the square of the average value of u[m]. 
Similarly, the total power in u[m] is equal to the average 
value of u 2 [m]. Thus, the power remaining for time- varying 
components of u[m] is: 

Avg(u 2 [m])-(Avg(u[m])) 2 =Avg((u[m]-Avg(u[m])) 2 ). 

This value is maximized by maximizing the dispersion of 
the samples about the mean. When the sample values are 
bounded, this maximization is achieved by placing half of 
the samples at each bound, so that the mean is equidistant 
from each bound. Since 0^u[m]^V4, this sets the mean of 
u[m] at Vs, and the maximum power present in harmonic 
components is: 

(.5*(0) 2 M).5*m 2 y(Vs) 2 =VM. 

Recall that at this worst case, half of the values of 
E{e 2 [m]} are zero and half are +!4. Since e 2 [m] is non- 
negative, E{e 2 [m]+0 implies that e[m]=0. Note that the 
difference between e[m] and e[m+l] is the phase increment 
modulo the quantization step size. If, for any m and m+1, 
e[m]=e[m+l]=0, the phase increment can be exactly 
expressed in the new quantization step. By induction, e[m] 
will be zero for all m if any two adjacent values E{e 2 [m]} 
and E{e 2 [m+1]} are both zero. The only possible sequence 
E{e 2 [m]} achieving the worst case is therefore 0, l A, 0, Va, 0, 
l A. . . This sequence has a single sinusoidal component at the 
Nyquist frequency, half the sampling rate. 

In the worst case the model to consider is u[m]=V&- 
(V^cosfam) since cos((27tf+7t)m+^)=cos((27tf-7t)m-Hl)). 
The expected waveform is: 

E{x[m}} - ((1 - te 2 A p 2/4) + 

(jc 2 A p 2/4)cos(7tm))cos(2jt^?i + O) + 0(A p 3) = 

(1 - it 2 A p 2/4)cos(27t>i + O) + 

(7C 2 A p 2/4)cos((27r/+ n)m + O) + 0(Ap3), 

clearly showing the desired signal and spur components. 
Thus, dropping the 0(A /7 3 ) term, a -18 dB per bit power 
behavior, the worst-case spur level relative to the desired 
signal after truncating to b bits is (SpSR): 


tc 4 A//16 JiTAn 

SpSR - -- -- o — 7 ^— - 7.84 - 12.042? dBc. 

(1 - n 2 A p 2 /4) 2 16 

This worst case is achieved for a large class of frequen- 
cies. Let the phase of the desired signal $=0, and let the 
desired frequency be R cycles in 2 &+1 samples, where R is 
an odd integer. The sequence e 2 [m] will be deterministic: 0, 
A p 2 /4, 0, Ap 2 /4. . . , exactly the worst case analyzed above. 
The frequency of the spur is the reflection of the desired 
signal across Va the sampling rate, and, as a result, it can be 
as close as Vi b th of the sampling rate from the desired 
sinusoid. 

In summary, if b bits of phase are output to a look-up 
table, and B bits of phase (B>b) are used prior to truncation, 
then the addition of an appropriate dithering signal using 
(B-b) bits will allow the word length reduction without 
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introducing spurs governed by the usual -6 b dBc behavior. 

If a single random variate is added as a dither signal 
(first-order dithering), the spur suppression is accelerated to 
12 dB per bit of phase representation. Since the table size is 
only effected linearly by the number of bits in a table entry, 5 
rather than exponentially as it is by the number of phase bits, 
the amplitude word length is of secondary importance to the 
phase word length, especially in all-digital systems. For 
example, -90 dBc spur performance would nominally 
require b=16 bits of phase and a 65,536 entry table. With 10 
first-order dithering, this level of performance requires only 
b 1 = 8.1 bits of phase per cycle in the look-up table address- 
ing. Worst-case spur performance of -100.5 dBc is achieved 
with 9 bits, a 512 entry table at most, and, at a 160 MHz 
sampling rate, TABLE I shows that with these realistic 15 
system parameters, the noise power spectral density is at a 
low -126 dBc/Hz. 

ACCELERATED SPUR SUPPRESSION 

20 

Further analysis based on an extension of results by Gray 
(R. M. Gray andT. G. Stockham, Jr., “Dithered Quantizers,” 
presented at 1991 IEEE International Symposium on Infor- 
mation Theory, Budapest, Hungary, June 1991) indicates 
that the phase spur suppression rate can be increased in steps 25 
of 6 dBc per bit by adding multiple uniform random deviates 
to the phase value prior to truncation. The addition of M 
uniform random deviates produces a dither signal with 
Mth-order zeroes in its characteristic function, thus making 
the Mth moment of the quantization error independent of the 3Q 
input sequence. 

An example of this technique providing 1 8 dBc per phase 
bit spur performance is shown in FIG. 6 . This technique 
involves adding two B-b bit uniform deviates to produce a 
B-b+1 bit dither signal, which achieves the accelerated spur 35 
reduction due to second-order zeroes in the dither charac- 
teristic function. Simulation results for when two uniform 
variates are added to the phase are presented in the next 
section. A straightforward extension of this technique to 
polynomial series allows spur-reduced synthesis of periodic 40 
digital signals with arbitrary waveforms. 

SIMULATION RESULTS 

Simulations were performed to validate the results of this 
analysis. These results were obtained using 8192-point 45 
un windowed FFTs, and the synthesized frequencies were 
chosen to represent worst-case amplitude and phase spur 
performance. FIG. 7 shows the power spectrum of a sine 
wave of one- eighth the sampling frequency truncated to 8 
bits of amplitude without dithering. " 

FIG. 8 shows the same spectrum with a sixteen-bit 
sinusoid amplitude dithered with one uniform variate prior 
to truncation to 8 bits. Note that the spurs have been 
eliminated to the levels consistent with those imposed by the 55 
initial sixteen bit quantization. 

FIG. 9 shows the spectrum of a 5-bit phase-truncated 
sinusoid with high-precision amplitude values. A worst- case 
example of first-order phase dithering is shown in FIG. 10. 
The measured noise power spectral density in FIG. 10 is 60 
-62.3 dBc per FFT bin, giving a noise density of -23.2-19 
log(fy/2) dBc, in agreement with the upper bound derived in 
Equation 7. The spur level is -52.3 dBc in the first-order 
dithered FIG. 10. 

FIG. 11 shows the same example using second-order 65 
(M=2) dithering using the sum of two uniform deviates. 
While the spectrum in FIG. 10 shows the residual spurs at 
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-12 dBc per bit due to second-order effects, FIG. 11 shows 
no visible spurs, indicating better than -63 dBc spurious 
performance. Additional simulations involving Megapoint 
FFTs and not represented by Figures confirm the -18 dBc 
per bit performance of the second-order phase-dithered 
system. 

Finally, FIG. 12 shows a worst-case result for first-order 
phase dithering together with first-order amplitude dithering. 
The amplitude samples are truncated to 8 bits, as are the 
phase samples. Note that the spurs are not visible in the 
spectrum; however, close analysis has demonstrated that 
they are present at the -88 dBc level expected due to 
second- order effects. 

A SYSTEM DESIGN EXAMPLE 

The block diagram of a direct digital frequency synthe- 
sizer based on the techniques presented here is shown in 
FIG. 13. The following system would perform at a sampling 
rate of 160 MHz, producing 8 -bit digital sinusoids spur-free 
to -90 dBc with better than -120 dBc/Hz noise power 
spectral density. The system parameters are as follows: 

Phase bits are in unsigned fractional cycle representation 
with: 

phase accumulator word-length determined by fre- 
quency resolution, and 

^16 bits prior to addition of 1 uniform phrase dither 
variate, with ^9 bits after dither addition and trun- 
cation; 

Amplitude look-up-table with: 

^2 7 =128 entries (using quadrant symmetries of ^16 
bits each normalized so that the sinusoid amplitude 
equals 512 16-bit quantization steps less than the 
full-scale value; 

Linear feedback shift register PN generator with ^16 lags 
producing one 8 -bit amplitude dither variate, and 

One LFSR PN generator with ^18 lags for generation of 
the 7-bit phase dither variate. 

CONCLUSION 

A digital dithering approach to spur reduction in the 
generation of digital sinusoids has been presented. A class of 
periodic dithering signals has been analyzed because of its 
similarity to LFSR PN generators. 

The advantage gained in amplitude dithering provides for 
spur performance at the original longer word length in an 
ideal system when the digital dithering signal is white noise 
distributed evenly, not uniformly, over one quantization 
interval. The reduced word length allows the use of fast, 
coarse-re solution, highly-linear digital-to- analog converters 
(DACs) to obtain sinusoids or other periodic waveforms 
whose spectral purity is limited by the D AC linearity, not its 
resolution. These results suggest that coarsely quantized, 
highly-linear techniques for digital-to -analog conversion 
such as delta- sigma modulation would be useful in direct 
digital frequency synthesis of analog waveforms. 

The advantage gained in the proposed method of phase 
dithering provides for an acceleration beyond the normal 6 
dB per bit spur reduction to a 6 (M+l) dB per bit spur 
reduction when the dithering signal consists of M uniform 
variates. Often the most convenient way to generate a 
periodic waveform is by table look-up with a phase index. 
Since the size of a look-up table is exponentially related to 
the number of phase bits, this can provide a dramatic 
reduction in the complexity of NCO’ s, frequency synthe- 
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sizers, and other periodic waveform generators. 

The advantages of dithering come at the expense of an 
increased noise content in the resulting waveform. However, 
the noise energy is spread throughout the sampling band- 
width. In high bandwidth applications, dithering imposes 5 
modest system degradation. It has been shown that high 
performance synthesizers with dramatically reduced com- 
plexity can be designed using the dithering method, without 
resulting in high noise power spectral density levels. 

Having thus described the invention in exemplary form, 10 
what is claimed is: 

1. An apparatus for generating a digital sinusoid signal 
having reduced spurious signal content; the apparatus com- 
prising: 

means for generating a digital sinusoid having a phase 15 
sample value and an amplitude sample value at an 
instant in time; 

a digital random noise signal generator; 
means for summing said noise signal with a parameter of 2 o 
said digital sinusoid to produce an N-bit word; and 
means for changing the wordlength of said N-bit word to 
produce an M-bit word where M is less than N; 
wherein said digital sinusoid parameter is said amplitude 
sample value of said sinusoid. 25 

2. An apparatus for generating a digital sinusoid signal 
having reduced spurious signal content; the apparatus com- 
prising: 

means for generating a digital sinusoid having a phase 
sample value and an amplitude sample value at an 30 
instant in time; 

a digital random noise signal generator; 
means for summing said noise signal with a parameter of 
said digital sinusoid to produce an N-bit word; and 35 
means for changing the wordlength of said N-bit word to 
produce an M-bit word where M is less than N; 
wherein said noise signal is summed with said phase 
sample value; said apparatus comprising a second 
digital random noise signal generator and means for 40 
summing said second noise signal with the amplitude 
sample value of said truncated M-bit word. 

3. The apparatus recited in claim 2 further comprising an 
amplitude look-up-table for generating said amplitude 
sample value corresponding to said phase sample value. 45 

4. The apparatus recited in claim 2 wherein said noise 
signal is the sum of K separate noise signals where K>1. 

5. An apparatus for generating a digital sinusoid signal 

having reduced spurious signal content; the apparatus com- 
prising: 50 

means for generating a digital sinusoid having a phase 
sample value and an amplitude sample value at an 
instant in time; 

a digital random noise signal generator; 5g 

means for summing said noise signal with a parameter of 
said digital sinusoid to produce an N-bit word; and 
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means for changing the wordlength of said N-bit word to 
produce an M-bit word where M is less than N; 

wherein said noise signal is the sum of K separate noise 
signals where K>1. 

6. A method for generating a digital sinusoid signal having 
reduced spurious signal content; the method comprising the 
steps of: 

a) generating a digital sinusoid having a phase sample 
value and an amplitude sample value at an instant in 
time; 

b) generating a digital random noise signal; 

c) summing said noise signal with a parameter of said 
sinusoid; and 

d) reducing the length of the sum resulting from step c); 

wherein step c) is performed by summing said noise 

signal with the amplitude value of said sinusoid. 

7. A method for generating a digital sinusoid signal having 
reduced spurious signal content; the method comprising the 
steps of: 

a) generating a digital sinusoid having a phase sample 
value and an amplitude sample value at an instant in 
time; 

b) generating a digital random noise signal; 

c) summing said noise signal with a parameter of said 
sinusoid; and 

d) reducing the length of the sum resulting from step c); 

wherein step c) is performed by summing said noise 

signal with the phase sample value of said sinusoid, the 
method further comprising the steps of generating a 
second digital random noise signal and summing said 
second noise signal with the amplitude sample value of 
said sinusoid. 

8. The method recited in claim 7 further comprising the 
steps of providing an amplitude look-up- table and generat- 
ing said amplitude sample value from said phase sample 
value using said look-up-table. 

9. The method recited in claim 7 wherein step b) com- 
prises the step of summing K separate noise signals where 
K>1. 

10. A method for generating a digital sinusoid signal 
having reduced spurious signal content; the method com- 
prising the steps of: 

a) generating a digital sinusoid having a phase sample 
value and an amplitude sample value at an instant in 
time; 

b) generating a digital random noise signal; 

c) summing said noise signal with a parameter of said 
sinusoid; and 

d) reducing the length of the sum resulting from step c); 

wherein step b) comprises the step of summing K separate 

noise signals where K>1. 
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